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The origin of small mixing among the quarks and a large mixing among the neutri¬ 
nos has been an open question in particle physics. In order to answer this question, 
we postulate general relations among the quarks and the leptonic mixing angles at 
a high scale, which could be the scale of Grand Unified Theories. The central idea 
of these relations is that the quark and the leptonic mixing angles can be unified at 
some high scale either due to some quark-lepton symmetry or some other underly¬ 
ing mechanism and as a consequence, the mixing angles of the leptonic sector are 
proportional to that of the quark sector. We investigate the phenomenology of the 
possible relations where the leptonic mixing angles are proportional to the quark 
mixing angles at the unification scale by taking into account the latest experimental 
constraints from the neutrino sector. These relations are able to explain the pattern 
of leptonic mixing at the low scale and thereby hint that these relations could be pos¬ 
sible signatures of a quark-lepton symmetry or some other underlying quark-lepton 
mixing unification mechanism at some high scale linked to Grand Unified Theories. 
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I. INTRODUCTION 


The quark mixing matrix, Vckm , parametrizes the misalignment in the diagonalisation 
of the up and down type quark mass matrices. It is well known that Vckm is almost close 
to a unit matrix. This implies that the quark mixing angles are small. On the other hand, 
analogous misalignment in the leptonic sector is encoded in the neutrino mixing matrix, 
Upmns ■ If turns out that Upmns is not close to a unit matrix. The mixing angles in the 
neutrino sector are large except 6 Vi [l-3j]. The origin of small mixing among quarks and a 
large mixing in the neutrino sector poses an intriguing open question. 

Among many approaches to explain the mixing pattern of the leptons, the assumption 
of family or flavor symmetries is a popular one. These symmetries differentiate among the 
members of different families and are usually discrete, finite and non-abelian, for reviews 


see Refs. 


,5]. 


This approach has been intensively used to study the mixing in the leptonic 


sector [6H9|. In addition to the leptonic mixing, there are also considerable efforts to un¬ 


derstand the quark mixing through family symmetries 


10 


The fami 


also be a built-in characteristic of the Grand Unified Theories (GUT) 13]. 


y symmetries can 


he quark-lepton unification is one of the most attractive features of the GUT theories 


14-1161]. The GUT symmetry group contains quarks and leptons in a joint representation. 


The weak interaction properties of the quarks and the leptons therefore get correlated. Hence 
it is possible in these theories, to derive the origin of the small and the large mixing in the 
quark and the lepton sectors respectively, along with any relation between them, if it exists. 

There are also reasons to speculate about the quark-lepton unification even on the exper¬ 
imental side. The so-called quark-lepton complementarity (QLC) relation jl7, 3 ] between 
the leptonic mixing angle 9 12 and the Cabibbo angle 9c 


7T 


( 1 ) 


012 + 0 C « 

can be a footprint of a high scale quark-lepton unification 
observation is due to the recent non-zero measurement of leptonic mixing angle # 13 22-1261 


17 


21]. Another interesting 
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which is 


9 


13 


9_c_ 

V2' 


( 2 ) 


This relation also hints a possible link between the quark and leptonic mixing, and it can 
be an artifact of some high scale quark-lepton symmetry in an underlying GUT theory [27]. 

Therefore, the present state of the measured leptonic mixing angles provide the theoretical 
motivation for a common origin of the quark and leptonic mixing at some high scale. In 
fact, the idea that the quark and lepton mixing can be unified at some high scale, referred 


to as “high scale mixing unification” (HSMU) hypothesis, was first proposed in Ref. 


28 


3l| . In recent studies 32-[34] it has been shown that HSMU hypothesis ‘naturally’ leads 
to nonzero and a small value for the leptonic mixing angle 0 13 and predicts a non-maximal 


623 ( cf. [32, [34| for details). This hypothesis has been studied in the context of Dirac 


neutrinos as well ( cf. 


35| for details). The central idea of this hypothesis is that the quark 

1 scale (referred to as the 
35] for details). In other 


mixing angles become identical to that of the leptons at some hig 


unification scale) which is typically taken as GUT scale (cf. 32 
words, at the unification scale 

O 12 = 0 \ 2 , d 13 = 9l 3 , 023 = 0f 3 , 


33 


(3) 


where 9ij (with i, j = 1, 2, 3) are leptonic mixing angles and 9\- are the quark mixing angles. 
This hypothesis nicely explains the pattern of mixing in the neutrino sector including the 


recent observation of nonzero and a small value of 0 13 


5-26] 


The large leptonic mixing 


angles at the low scale are obtained through the renormalization group (RG) evolution of 
the corresponding mixing parameters from the unification scale to the low scale. 

The implementation of the HSMU hypothesis requires the minimum supersymmetric 
standard model (MSSM) as an extension of the standard model (SM). The working of the 
HSMU hypothesis is as follows. We first evolve the quark mixing angles from the low scale 
(mass of the Z boson) to the supersymmetry (SUSY) breaking scale using the SM RG 
equations. After that, from the SUSY breaking scale to the unification scale, evolution of 
quark mixing angles is governed by the MSSM RG equations. I 11 the next step, the quark 
mixing angles at the unification scale, are put equal to that of the neutrinos following the 
HSMU hypothesis. The leptonic mixing parameters are then run from the unification scale 
to the SUSY breaking scale using the MSSM RG equations. From the SUSY breaking scale 
to the low scale, mixing parameters are evolved through the SM RG equations. 




























4 


In addition to SUSY, we also need a large tan/3 to realize the HSMU hypothesis [32|, 33]. 


The only free parameters during the top-down running of the leptonic mixing parameters 
are masses of the three light neutrinos. They are chosen at the unification scale in such a 
manner that we recover all the mixing parameters at the low scale within the 3cr limit of the 
global fit. It turns out that the chosen masses of neutrinos must be quasi-degenerate (QD) 


and normal hierarchical 


32 


33] 


In this work, inspired by the HSMU hypothesis, we postulate the most general relations 
among the quark and the leptonic mixing angles at the unification scale. In a compactihed 
form the most general relation among the leptonic and the quark mixing angles within the 
same generations is as following 


$i2 = a? 1 $? 2 , 9 13 = a k 2 2 9l 3 , 023 = 0^3- (4) 

where k %, with i = (1,2,3) are real exponents. We refer to this relation as the “high scale 
mixing relation’’ (HSMR). We have chosen (/ci, fc 2 , £ 3 ) to be (1,1,1) for the simplicity of our 
analysis. The relations within the same generations are the simplest generalization of the 
HSMU hypothesis. In principle, we can also construct the most general HSMR relations 
among different generations completely independent of the HSMU hypothesis. The analysis 
of these relations is beyond the scope of this work and could be studied elsewhere. 

There will be different possibilities depending on the relations among the proportionality 
factors. We firstly list below the different possible cases with the maximum and the minimum 
allowed values of the three independent proportionality factors cq, 


Case A : 

$12 

= < ax $? 2 , 

$13 

— n/ max O q 

— CK 2 17 13) 

$23 

= «r Di>, 

(5) 

Case B : 

$12 

= $?2) 

$13 

_ man nq 

~ a 2 $13> 

$23 

= af" <4. 

(6) 

Case C : 

#12 

_ max nq 
~ <*l #12’ 

$13 

_ „,min nq 

~ a 2 $13) 

$23 

= a?”* 023. 

(7) 

Case D : 

$12 

— ry max f) q 

— Uf 1 (7 12 , 

$13 

— cy min 
~ u 2 u lZ-> 

$23 

= <" 0 I 3 . 

(8) 

Case E : 

$12 

_ min nq 
= <*1 #12’ 

$13 

„max nq 
= «2 $13) 

$23 

~ a 3 #23’ 

(9) 

Case F : 

$12 

= of" 0? 2 , 

$13 

= «2 $13) 

$23 

= < n 0f 3 . 

(10) 

Case G : 

$12 

= < ” 0?2. 

$13 

— n/ min f) q 
~ a 2 *13> 

$23 

= ar* 023. 

(11) 

Case H : 

$12 

= of” e? 2 . 

$13 

_ „min nq 
~ a 2 $13) 

$23 

= «r 

(12) 


In this work, we have presented our results for the maximum and minimum allowed values 
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of a, for all the above cases, Eqs. (l51 fT2|) . We then move on to scenarios, assuming relations 
among the cq’s. There can be more general HSMR where two proportionality constants 
can be identical and the third one is different. However we will discuss in this work more 
simplified scenerios, where the three proportionality constants are equal. 

$i 2 = « fcl $? 2 , e 13 = a k2 9l 3 , d 23 = c/ 3 $ 2 V (13) 

As explained before we have restricted to values of hi as either 0 or 1 . We note that the value 
(ki, k 2 , k 3 ) = (0, 0, 0) will reduce HSMR to HSMU hypothesis making Eq. ((3]) a specific form 
of HSMR, Eq. (TT3|) . We present below the seven different possible cases, where the quark 
mixing angles are assumed to be proportional to the corresponding leptonic mixing angles. 


Case 

1 : 

9l2 

= « $12 > 

$13 

= $13 > 

$23 

CO 

(14) 

Case 

2 : 

$12 

~ $12) 

$13 

= a $13) 

$23 

= $23) 

(15) 

Case 

3 : 

$12 

= $12) 

$13 

= $13 > 

$23 

= a $23) 

(16) 

Case 

4 : 

$12 

= a $i 2 > 

$13 

= a 6f 3 , 

$23 

= $23) 

(17) 

Case 

5 : 

$12 

~ $12) 

$13 

= a $ 13 ) 

$23 

= a $23) 

(18) 

Case 

6 : 

$12 

= a $12 j 

$13 

= $13) 

$23 

= a $23) 

(19) 

Case 

7 : 

$12 

= a $ 12 ) 

$13 

= « $i 3 ) 

$23 

= « $23- 

(20) 


The proportionality constant a in the above Eqs. (I14H20P is taken as real parameter. We 
have carried out a detailed study for these cases in this work. 

We note that there exist GUT models in the literature where proportionality between the 
quark and the leptonic mixing angles are explicitly shown. For example, the proportionality 
relation observed between the leptonic mixing angle $13 and the Cabibbo angle 9c in Eq. (|2]) 
can arise naturally in SU(5) GUTs and Pati- Salam models. For more details, see [27]. 


Further more, it is shown in Ref. 


18] that the relations between the quark and the leptonic 


mixing angles are possible and they support the idea of grand unification. However, non- 
abclian and abelian flavor symmetries are essential to make this happen 18] . 

There is two-fold importance of the HSMR hypothesis. The first remarkable feature 
is that these relations provide a very simple way to achieve a large neutrino mixing. We 
shall see that predictions of these relations are easily testable in present and forthcoming 
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experiments. The second importance is that if predictions of HSMR hypothesis are confirmed 
by experimets, like neutrinoless double beta decay, this would be a strong hint of quark- 
lepton unification at high scale. 

The plan of the paper is as follows. In section [Til we present the required RG equations 
for the running of the neutrino mixing parameters. The SUSY threshold corrections and 
the neutrino mass scale are discussed in section [1111 The results are presented in section [TV] 
using dimensional-5 operator as well as in the framework of type-I seasaw. In section [Vj for 
the sake of illustration, we discuss two models where HSMR hypothesis can be realised. We 
summarize our results and conclude in section ED 


II. RG EVOLUTION OF THE LEPTONIC MIXING PARAMETERS 


In this section, we briefly discuss the RG evolution of the leptonic mixing parameters.The 
most often studied scenario is the one where the Majorana mass term for the left handed 


neutrinos is given by the lowest dimensional operator 


36] 


X = jKj/CVV.< + h ' C ' • 


( 21 ) 


in the SM. In the MSSM, it is given by 


^ISSM = W J + = _i L^hfl + h.c. , 


( 22 ) 


where k 9 j has mass dimension —1, tf;' is the charge conjugate of a lepton doublet and 
a,b,c,d G {1,2} are SU(2) L indices. The double-stroke letters 1L and lh denote the lepton 
doublets and the up-type Higgs superfield in the MSSM. Using this mass operator, we 
introduce neutrino masses in a rather model independent way since it does not depend on 
the underlying mass mechanism. 


The evolution of the above dimensional-5 operator below the scale w’ 


is provided by its RG equation. The one loop equation is as follows 


16t r 2 k = C (YjY e ) r k + Ck (YjY e ) +&k, 


37-4' 


rere it is generated 


( 23 ) 


where k = t — ln(/i//xo) and /i is the renormalization scale and 


C = 1 


in the MSSM , 





in the SM . 
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C =-2 


(24) 


The parameter a in the SM and MSSM is given by 

«sm = — 3#2 + 2 (y~ + y\ + He) + 6 (?/* + 2/b + 2/c + 2/s + 2/d + 2/«) + ^ > 
6 

«mssm = — -Pi — 6^2 + 6 (y f 2 + y'c + y'y) ■ 


(25) 


The quantities yj (/ G {e, d, u}) represent the Yukawa coupling matrices of the charged 
leptons, down- and up-type quarks respectively, g t (i = R 2) denote the gauge couplings and 


A is the Higgs self coupling. For more details see Ref. 


36]. 


We are interested in the RG evolution of parameters that are the masses, the mixing 
angles and the physical phases. The mixing angles and the physical phases are described by 
the PMNS matrix. This matrix is parameterized as follows 


U PM NS — V • U, 


(26) 


where 


/ 


V = 


C 12 C 13 


S 12 C 13 


si 3 e 


—i<5 \ 


-C23S12 — 5235l3Cl2e 1<5 C23C12 — S23’Sl3' S 12e 1<5 S23C13 


,i<5 


i<5 


S23S12 — C23Sl3Ci2e —S23C12 — C23Si3Si2e C23C13 


(27) 


and 


/ e -bu/ 2 


U = 


V 


0 

0 


0 0 

e -iw/2 0 

0 1 


\ 


with Cij and s tJ dehned as cos % and sin% (i,j = 1,2,3), respectively. The quantity S 
is the Dirac phase and (p±, ip 2 are the Majorana phases. The global experimental status of 
the leptonic mixing parameter is summarized in Table Q] 

Here we would like to remark that the RG equations (I23p for Yukawa couplings matrices 
are parametrization independent. The main aim is to probe if there is any connection 
between the quark and the leptonic mixing. For this purpose, we have chosen the standard 
parametrization which is the most studied and also commonly used in the literature. In 
principle, one could use an alternative parameterization to work and test the reality of 
HSMR. The results can be always interpreted as a possible indication of a connection between 
quark and leptonic mixing. 







Quantity 

Best Fit 

3cr Range 

(10 -5 eV 2 ) 

7.54 

6.99 - 8.18 

A m 32 (10 -3 eV 2 ) 

2.39 

2.20 - 2.57 

0?2 

33.71 

30.59 - 36.81 

023 

41.38 

37.7 - 52.3 

013 

8.8 

7.63 - 9.89 


TABLE I: The global fits for the neutrino mixing parameters 


Q. 


We now summarize the RG equations used for running the leptonic mixing par ameters 
from high to the low scale. For a detailed discussion of these equations, see Ref. [36|. These 
equations are derived using the lowest dimensional neutrino mass operator as discussed 
above and are given by the following analytical expressions 36] 


012 — 


Cyl . nn 2 + m 2 e ilf2 \ 2 

sin 2 ui 2 ^23 -i—o- 


32tt 2 


Ami 


131 


21 


(28) 


Cv 2 

^13 = sin 2012 sin 2023 


m 3 


x 


32vr 2 - — Am| 2 (l + C) 

x [mi cos((pi - 5) - (1 + () m 2 cos(^ 2 - S) - (m 3 cos 5] + ^(0i 3 ) , (29) 


023 — 


Cy * sin 202 3 1 


32?r 2 


Am§ 2 L 


q 2 |m 2 e 1 ^ 2 + m 3 |- + s 


2 2 |mi e lipi + m 3 | 2 


12 


1 + C 


13 ) 


(30) 

where 0y = ^ (with i , j = 1, 2, 3), t = ln(p//z 0 ), /i being the renormalization scale and 

(31) 


C : = 


A ml 


_21 

Am § 2 


, A/7% := m^ — mf, A m 32 := mi — m^. 


For the masses, the results for y e = y Pb = 0 but arbitrary 0 33 are 


167r 2 mi = 

[a + Cy 2 r ( 2 s 2 12 s 23 + F\)] m 3 , 

(32a) 

167r 2 m 2 = 

[a + Cy 2 T ( 2 c 2 2 s\ 3 + F 2 )] m 2 , 

(32b) 

167r 2 m 3 = 

[a + 2Cy 2 cl 3 c 2 3 ] m 3 , 

(32c) 


where m, = (i =1, 2, 3) and Fi, F 2 contain terms proportional to sin0i 3 , 

F x = — s i3 sin 2012 sin 20 23 cos S + 2s 2 3 c 2 2 (? 23 , 


(33a) 
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F 2 = si 3 sin 2012 sin 2023 cos $ + 2s ^ 3 s ^ 2 C 23 . (33b) 

In this work, we are working in the CP conserving limit which means Majorana and Dirac 
phases are assumed to be zero. Therefore, we have not provided the RG equations for 
them. The non-zero phases are expected to have non-trivial impact on the parameter space. 
However, this study is beyond the scope of the present work and will be presented in a 
future investigation. Furthermore, we also study the effect of the new physics which could 
generate the above dimensional-5 operator. For this purpose, we present our analysis within 
the framework of type-I seesaw. 

Now, we briefly discuss the evolution of the leptonic mixing angles. In the SM as can be 
seen from Eq. (1251) , only tau Yukawa coupling will dominate the evolution which is already 
very small. Hence the running of the neutrino masses is governed by a common scaling factor 
and the evolution of leptonic mixing angles can only be enhanced for QD mass pattern. In 
the MSSM the value of tau Yukawa coupling can be larger with respect to the value in the 
SM for a large value of tan /3. Hence the evolution of the leptonic mixing parameters can 
be enhanced in addition to the enhancement coming from the QD neutrino mass pattern as 
discussed below. 

It is interesting to note from Eqs. (1251 [251 and 1501) that the major contribution to RG 
evolution of the mixing angles arises due to following enhancement factors 

#12 oc £ 1 , 0i3, 0 23 oc (34) 

where 

t _ m2 t _ m2 /or\ 

1 A m 2 21 ’ 2 Am § 2 ’ 

and m is the average neutrino mass with m = (rrii + m 2 + m^)/2>. It is clear that we need 
masses of the neutrinos to be QD to explain the largeness of mixing angles at the low scale. 

III. THE LOW ENERGY SUSY THRESHOLD CORRECTIONS AND THE 
ABSOLUTE NEUTRINO MASS SCALE 

We discuss the required low energy SUSY threshold corrections for the mass square 
differences and the significance of the absolute neutrino mass scale in this section. 




10 


A. The low energy SUSY threshold corrections 


It is well established in the previous works on HSMU hypothesis that among the five 
mixin g pa rameters, one of the mass square differences (Am^) lies outside the 3 cr global 


range 


28 


32], As shown in the previous works, this mass square difference can be brought 


well within the 3cr global limit, i: 
to the mass square differences 


the low energy SUSY threshold corrections are incorporated 


28 


33. The importance of SUSY threshold corrections for 


QD neutrinos is discussed in Refs. 41n44]. These corrections are given by the following 


equations 


29] 


(A m 21 )th 

l 32 )th 
Hl)th 


(A 177.32 
(A 777.3 


2m 2 cos 2 ^ 12 ]—2T e + T^ + T r ], 
2m 2 sin" 0 12 [—2T e + T M + T r ], 
2m 2 cos 2 612 [—2T e + T^ + T r ]. 


(36) 


given by 


is 

le 

41, 

44 


92 


32vr 2 [ y^Va 


x l “ 4 + (# 


i) 


y& 


ln(x\ 


^ ^ ln(xD 

vl 


(37) 


where g 2 is the 577(2) coupling constant and y& — 1 — x\ with = M&/M.M.$ stands 
for wino mass, M& represents the mass of charged sleptons. We work with an inverted 
hierarchy in the charged-slepton sector where the mass of selectron is defined through the 


ratio R = tt^-. The mass of the wino is chosen to be 400 GeV following the direct searches 


at the LHC 


45]. 


B. The absolute neutrino mass scale 

The scale of the neutrino mass is one of the open questions, ever since it has been 
confirmed that the neutrinos are massive. I 11 case of QD and the normal hierarchical spectra, 
we have 

mi<m 2 <m 3 ~ m 0 (38) 

mo \J A 77732 ~ 5 x 10~ 2 eV. 


with 


(39) 
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There are three complementary ways to measure the neutrino mass scale. The first one, a 
model independent method, is to use the kinematics of /3-decay to determine the effective 
electron (anti) neutrino mass imp). It is given by 


771/3 = 


£i^ 


I 2 777 2 


(40) 


The 777/3 has an upper bound of 2 eV from tritium beta decay 


46 


47|. In future, the KATRIN 
experiment has sensitivity to probe mp as low as 0.2 eV at 90% CL [48]]. We note that m 0 in 
the QD regime for CP conservation is approximately equal to the effective beta decay mass 
mp. Hence QD mass pattern is well within the sensitivity of the KATRIN. 

The second method to extract the neutrino mass is neutrinoless double beta decay which 
assumes that neutrinos are Majorana particles 49j, [50]. The observable parameter M ee , the 


double beta decay effective mass is given as following 


M ee = \^U 2 ei mi 

= 177i 1 cf 2 c 2 3 e~ llpl + 777 2 s^ 2 c 2 3 e ~ llp2 + m 3 s\ 3 e~ l25 \ 


(41) 


For quasi-degenerate neutrinos 


M ee ~ 777 0 |c 2 2 c 2 3 e i(pi + s 32 c 33 e llfi2 + s 33 e 


2 -i25 1 


(42) 


Since the contribution of m 3 is suppressed by the small sin 2 0 13 coefficient, we obtain 


M ee ~ 777,0 \/1 — sin 2 2$i 2 


2 07 ) (1 - cos(^i - </? 2 )) 


(43) 


For CP conserving case where the Majorana and Dirac phases are zero, M ee ~ mo- For 
M ee ~ 0.1 eV, the above expression corresponds approximately to half-life in the range 
of 10 25 to 10 26 yrs 49| which makes the QD mass scheme testable in present and future 


experiments . In the QD regime, the neutrino mass can be written as 


49] 


mo < (■ M ee )l ^ 


1 + tan 2 612 


= (M ee )^ x f( 6 12 : 6 13 ). 


(44) 


/max 1-tan 2 e u -2\U e3 \ 2 
Using inputs from Table [U the function /(#i 2 ,#i 3 ) has a range from 2.2 to 4.1 at 3cr. The 
most stringent upper limit on the effective mass M ee provided by the GERDA experiment 


is 0.4 eV 


511 ]. Hence m 0 < 1.64 eV and sum of the neutrino masses £ 777 , = 3r77 0 < 4.91 eV. 


















12 


The third determination of neutrino masses is provided by the cosmological and astro- 
physical observations. The sum of the neutrino masses, Em,, has a range for upper bound 
to be 0.17 — 0.72 eV at 95% CL 52], This limit is not model independent and depends on 


the cosmological model applied to the data. 


IV. RESULTS 


We present our results in this section for the different cases listed in Eqs. (in - EDI) 
and for limiting cases of the most general HSMR as shown in Eqs. As discussed 

earlier, we need MSSM as an extension of the SM for the implementation of HSMR and 
HSMU hypothesis. In the first step, we run quark mixing angles, gauge couplings, Yukawa 
couplings of quarks and charged leptons from the low scale to the SUSY breaking scale. The 
evolution from the SUSY breaking scale to the unification scale is done through the MSSM 
RG equations. After evolving up to the unification scale, we obtain quark mixing angles 
6\ 2 = 13.02°, 0f 3 = 0.17° and 6 ^ = 2.03°. In the next step, quark mixing angles are used 
to calculate the leptonic mixing angles using HSMR at the unification scale. After this, we 
run down the MSSM RG equations up to the SUSY breaking scale. The SM RG equations 
take over the evolution of mixing parameters beyond the SUSY breakin g sc ale. The SUSY 
breaking scale is chosen to be 2 TeV following the direct LHC searches [45|. We also need 
a large tan [3 which is chosen to be 55. The unification scale where HSMR can exist is 


chosen to be 10 14 GeV which is consistent with present exp erimental observations jl|. We 


have used the MATHEMATICA based package REAP 53j] for the numerical computation 
of our results. We have done a rigorous, thorough and comprehensive in this work. For this 
pupose, we have written an interface code which together with public code can be used to 
scan whole parameter space. 


A. RG evolution of HSMR 

We study the RG evolution of HSMR as given in Eqs. (1141 - 1201) and compare our results 
with respect to the HSMU hypothesis. In Fig. [T] we show how enhancement factors and 
£2 evolve from the unification scale to the low scale, as a deviates from unity. The results 
are displayed for all the cases. It can be seen from Fig. [T| that the evolution in Case 1 at 
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HSMU hypothesis 


Case 1 


Case 2 




- 5 , 


10000 1e+08 le-t-12 
(i (GeV) 



Case 3 


Case 4 


Case 5 




FIG. 1: The change in the RG evolution of the enhancement factors (£ 1 , 2)1 ESI for the different 
cases of HSMR as a function of the RG scale p when a deviates from unity. 

a = 1.1 is similar to HSMU hypothesis. However, as a approaches to lowest value on the 
left panel of Case 1, £1 changes sufficiently. Similarly for the upper limit of a = 1.962, the 
evolution again becomes very different from the HSMU hypothesis. This explains why the 
RG evolution of the PMNS mixing angles change when a deviates from unity. The same 
argument follows for all the other cases of HSMR and can be checked from Fig. Q] 

We next show the evolution of the mixing angles for the different cases in Fig. [2] along 
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HSMU hypothesis 


Case 1 


Case 2 



n(GeV) 




Case 3 


Case 4 


Case 5 





|i (GeV) |i(GeV) n(GeV) 


Case 6 Case 7 



FIG. 2: The RG evolution of the PMNS mixing angles with respect to the RG scale p, for the 
different cases, with a varied in the respective allowed range. 


with the HSMU hypothesis. We observe from Figs. |T| and [21 that the evolution of HSMR is 
similar to the HSMU hypothesis when a deviates slightly from unity. However, when a is 
very far from unity, RG evolution undergoes dramatic changes. There is another interesting 
phenomenon that can be observed from Fig. [2j It can be easily seen that the RG evolution 
of the mixing angles, for Cases 3 and 5 are similar, with #12 and $23 almost similar at the 
low scale at the lower end of a. The difference between them at the low scale increases with 
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the increase in value of a. The pattern is exactly opposite in the other cases of HSMR, with 
the difference between 9 12 and 6^3 at the low scale decreasing as one goes from the lower to 
the upper end of a. This in a way tells us beforehand that the phenomenology of Case 3 
and 5 will be similar, which will be discussed in detail afterwards. 

B. Phenomenology of HSMR 

In this subsection, we discuss in details the phenomenological implications of HSMR. 
Our aim is to investigate the behavior of a as it deviates from unity and its phenomeno¬ 
logical consequences taking into account all the experimental constraints of Table Q ] and the 
GERDA limit 5l|. The common observation among all HSMR is the emergence of the 


strong correlations among AmL M ee , 6 * 23 , 0 13 and Em* 


1. HSMU hypothesis 


As observed earlier, the value a — l will reduce all cases of HSMR to HSMU hypothesis. 
We present a full parameter scan of the HSMU hypothesis using dimensional-5 operator. It 


shoulc 


32 


be noted that this analysis was absent in the previous works on HSMU hypothesis 


34] and is reported in this work for the first time. We present a correlation in Fig. [3l 
which is not studied in the previous investigations. We show here the variation of A m \ 2 
with respect to M ee . The M ee has an upper bound of 0.4 eV from the GERDA experiment 
5ll ]. Using this limit, we are able to put an upper bound on the allowed range of A m\ 2 - The 


allowed range for A m \ 2 is (2.21 — 2.45) x 1CU 3 eV 2 as observed from Fig. [31 The lower bound 
on M ee is 0.384 eV for the HSMU hypothesis. Hence, our work on the HSMU hypothesis 
will be ruled out if GERDA crosses this number in the future. The effective /3 decay mass 
mp is another interesting observable since it does not depend on whether the neutrinos are 
Majorana or Dirac. The prediction for mp coincides with the effective double beta decay 
mass M ee in the QD regime and for CP conservation. Hence, the allowed range for mp is 
identical to that of M ee in our work. 

In Fig. [4j we show the variation of 623 , with respect to 613 . We observe a strong correlation 
between 0 13 and $ 23 - The difference between this investigation and that of presented in 


Ref. 


32] is the variation of 612 . In the previous work, this correlation was reported for a 
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FIG. 3: The variation of Am| 2 with respect 
to M ee , in the context of the HSMU hypoth¬ 
esis, with a = 1. 



FIG. 4: The variation of 0 2 3 with respect to 
013, in the context of the HSMU hypothesis, 
with a =1. 



chosen value of the angle 0 12 at the low scale in the context of type I seesaw. In this work, 
we do not choose any particular value of 0 12 at the low scale. We obtain a band for this 
correlation and previous results are a specific case of our present results. We observe that 
023 is non maximal and always lies in the second octant. This confirms the predictions of 


our earlier work 


32], The allowed range of 0 13 is 7.63° — 8.34° and that of 0 2 3 is 49° — 52.3 C 



M ee (eV) 


FIG. 5: The variation of 0i 2 with respect to 
M ee , in the context of the HSMU hypothesis, 
with a =1. 



FIG. 6: The variation of Era; with respect to 
M ee , in the context of the HSMU hypothesis, 
with a =1. 


We next present variation of 0 12 against M ee in Fig. 0 This correlation is also a new 
prediction of our work and do not exist in previous studies. The whole 3cr global range for 
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the angle #12 is allowed for the M ee < 0.4 eV, However, as can be observed from Fig. O 
the range 34.4° < $ 12 < 36.81° is ruled out for 0.384 eV < M ee < 0.393 eV. The precise 
predictions for all observables are provided in Table Hill In the end, we also have a new 
correlation between the sum of neutrino masses and M ee which is not studied previously. 
This correlation is shown in Fig. |6j Our prediction for sum of neutrino masses is 1.16 — 1.2 
eV using the upper bound on M ee given by the GERDA. 

2. The most general HSMR within the same generations 

The most general HSMR within the same generations for (/c 1; k 2 , Ay) = (1,1,1) as defined 
before is given by the following equation 

6*12 = on d\ 2 i $13 = «2 $ 13 , $23 = $| 3 . (45) 

We present the results for the maximum and the minimum values of a* for Eqs. fED, 
taking into account all the experimental constriants. In the Table [Til we present the allowed 
values of a* along with the the respective physical masses and the mixing angles. 

It is remarkable that in the Case E, all the mixing parameters are within 3cr global range 
without adding threshold corrections. If we add threshold corrections, the predictions are 
A7 ttJ 2 = 2.35 x (10 - 3 eV 2 ) and 7.01 x (10 _ 5 eV 2 ) for R = 1.0. Thus, threshold 

corrections at this point are effectively negligible. We further notice that the different 
combinations of the allowed end points of cq, leads to M ee around 0.35 eV - 0.4 eV. This 
most general case with different a*, alone will not suffice, when the value of M ee will be 
further constrained by the future experiments. We then have to look for more specific cases, 
where the ads will not be different, but have some relations among them. We consider the 
simplified scenario, where the ads are equal. We have carried out a detailed analysis for all 
the possible seven cases in this scenario in the next subsections. 

3. Case 1: 0\ 2 = a 0\ 2 , $13 = $135 $23 = $23 

The first case of HSMR is the one where leptonic mixing angle $12 is proportional to 
$^ 2 and the other two angles are identical. In Fig. [3 we show how the correlation between 
Am 2 2 and M ee changes as a deviates from unity. We observe on the left panel of Fig. [7] 
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Qtl 

«2 

<*3 

Masses at unification scale (eV) 

£m.j (eV) 

^12 

^13 

^23 

A m 2 32 (l<T 3 eV 2 ) 

A ml, (10~ 5 eV 2 ) 

M„ (eV) 

Lightest neutrino mass: mi (eV) 

R, 

mi 

m 2 

m 3 

Case A 

1.46 

2.54 

1.19 

0.4583 

0.461 

0.5186 

1.16 

36.52 

9.88 

41.14 

2.5 

8.06 

0.385 

0.3850 

2.29 

Case B 

1.45 

1.68 

0.91 

0.4757 

0.478 

0.5380 

1.20 

30.61 

8.79 

37.97 

2.25 

8.12 

0.40 

0.3997 

1.8 

Case C 

1.38 

0.71 

1.28 

0.489 

0.493 

0.5527 

1.24 

36.8 

9.87 

50.83 

2.22 

8.14 

0.411 

0.4106 

5.3 

Case D 

1.14 

0.92 

0.94 

0.4754 

0.478 

0.5370 

1.20 

31.18 

7.70 

45.31 

2.20 

8.14 

0.40 

0.3994 

1.69 

Case E 

0.8 

2.2 

1.15 

0.4096 

0.412 

0.4625 

1.04 

32.77 

7.65 

48.13 

2.35 

7.01 

0.344 

0.3442 

- 

Case F 

0.89 

1.61 

0.82 

0.4751 

0.477 

0.5361 

1.20 

30.6 

7.65 

43.66 

2.22 

7.37 

0.40 

0.3993 

1.06 

Case G 

0.92 

0.98 

1.03 

0.4421 

0.445 

0.4989 

1.12 

32.37 

7.64 

52.19 

2.22 

7.86 

0.372 

0.3714 

1.48 

Case H 

0.88 

0.95 

0.86 

0.4764 

0.479 

0.5372 

1.20 

30.99 

7.63 

51.97 

2.22 

7.55 

0.40 

0.4003 

1.29 


TABLE II: The allowed predictions for the different cases of the most general HSMR for minimum 
and maximum allowed values of cu, Eqs. mm- 

that the lowest allowed value of a is 0.902. This value is derived by the 3cr global limit of 
the leptonic mixing angles. On the right panel of Fig. [2 the upper bcmnd on a is shown. 
For the upper bound on a , in principle, one can go up to 1.962 with all mixing parameters 
within the global range. This value of a belongs to M ee > 0.4 eV and hence is ruled out 
by the GERDA limit. The allowed upper bound on a is 1.28 which is derived using the 
GERDA limit. 

We compare Fig. [7]with Fig. [B] of the HSMU hypothesis (a = 1) to study the phenomeno¬ 
logical behavior of a. As obvious from the left panel of Fig. [7]. M ee has its maximum allowed 
range at the lowest value of a. This is because the absolute neutrino mass decreases for 
a < 1 and increases for a > 1 in the case under study. Hence, at a = 0.902 on the 
left panel of Fig. [TJ we obtain 0.365 eV < M ee < 0.40 eV corresponding to whole 3<r 
global range of AThe same prediction for the HSMU hypothesis in Fig. [31 (« = 1) 
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FIG. 7: The variation of Am| 2 with respect 
to M ee for Case 1 of HSMR. 


FIG. 8: The variation of 023 with respect to 
013 for Case 1 of HSMR. 


is 0.385 eV < M ee < 0.418 eV which belongs to A m\ 2 = (2.21 — 2.45) x 10 -3 eV 2 . The 
prediction when a slightly deviates from unity (a = 1.1) is 0.384 eV < M ee < 0.435 eV 
corresponding to AmJ 2 = (2.22 — 2.57) x 10” 3 eV 2 . At the upper allowed value of a = 1.28, 
we have 0.4 eV < M ee < 0.45 eV which belongs to A m\ 2 = (2.20 — 2.57) x 10” 3 eV 2 . We 
observe that the uppermost value of a = 1.962 has 0.571 eV < M ee < 0.625 eV belonging 
to Am 2 2 = (2.23 — 2.57) x 10” 3 eV 2 . This value of a is already ruled out by the GERDA 
limit. 

This case can be ruled out if GERDA reaches M ee < 0.365 eV. There is an apparent 
overlap between predictions of the case under study and the HSMU hypothesis. This can 
be discriminated using the SUSY ratio R. For a clear picture of the phenomenological 
consequences, we provide values of mixing parameters and other observables belonging to 
minimum and maximum allowed values of a for each case and the HSMU hypothesis in 
Table ED 

The variation of 023 with respect to $13 is shown in Fig. [ 8 j The mixing angles reach their 
3cr limits at their lower and upper ends. For example, at a = 0.902, 0 13 is at its minimum of 
the 3 <t global limit while 0 23 is at its maximum independent of the upper bound of M ee . On 
the other hand, at a = 1.962, the predictions are reversed but this value is already rejected 
by the GERDA limit of M ee . The allowed ranges of 6*13 and 023, at a = 1.1, are 7.62° — 9.1° 
and 46.09° — 52.2° respectively. Compared to Fig. [4] for the HSMU hypothesis where the 
allowed range of 623 is always in the second octant, 023 has its minimum value 44.04° at 
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M„(eV) 


FIG. 9: The variation of Q\ 2 with respect 
to M ee for Case 1 of HSMR. 


FIG. 10: The variation of Em, with respect 
to M ee for Case 1 of HSMR. 


a = 1.28 which belongs to M ee = 0.4 eV and lies in the first octant. The corresponding 
value of 0 i 3 is 8.16°. 

We next show the behavior of # 12 with respect to M ee in Fig. [HI We observe that at 
a = 0.902, 612 is at its global minimum 30.6°. On the left panel of Fig. [HI at a = 1.1, 0 12 has 
an allowed range of 30.6° — 35.65° for M ee < 0.4 eV. For the HSMU hypothesis in Fig. 0 0 i2 
has the whole 3<r global range with some higher values ruled out for M ee < 0.393 eV. For 
a = 1.28, the value of d\ 2 is 32.82° for M ee = 0.4 eV as can be seen from the right panel of 
the figure. For a = 1.962, 0 12 reaches to the maximum of its 3cr global limit. 

Finally in Fig. [TUI the variation of sum of the neutrino masses with respect to M ee is 
presented. For lowest value of a = 0.902, on the left panel, the range of Em; is 1.12 — 1.2 eV 
for M ee < 0.4 eV. At a = 1.1, it is 1.16 — 1.2 eV for M ee < 0.4 eV. On the right panel, the 
value of sum at a = 1.28 is 1.2 eV and the region Emj >1.2 eV belongs to M ee > 0.4 eV. 


4■ Case 2: Q\ 2 = #13 = « #23 = #23 

The second case which we consider has leptonic mixing angle #13 proportional to In 
this case, the lower bound on a is 0.45 which is derived using global limits on the mixing 
angles. The a on the upper side, however, is remarkably bounded by the ratio R. This 
theoretical bound arises because we work with an inverted hierarchy in the charged-slepton 
sector and at a = 2.5, we have R = 1. In principle, a has a range up to 3.5 satisfying all 
experimental constraints with R < 1. 
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In Fig. dH we show the behavior of A m\ 2 versus M ee for different values of a. For the 
a = 0.45 on the right panel, we have 0.382 eV < M ee < 0.418 eV which corresponds to 
the whole 3cr global range of Am§ 2 . In case of a = 1.1 on the left panel, 0.38 eV < M ee < 
0.428 eV belongs to A m\ 2 = (2.2 — 2.57) x 10” 3 eV 2 . The range of M ee at the upper end 
a = 2.5 is 0.342 eV < M ee < 0.378 eV corresponding to Am 2 2 = (2.2 — 2.53) x 10 -3 eV 2 . 
A remarkable feature emerges in this case. Unlike case 1, the absolute neutrino mass scale 
increases for a < 1 and decreases for a > 1. Now, at the upper allowed value of a = 2.5, 
M ee is sufficiently below the GERDA limit. We would emphasize that one of the main 
observations of this case is that a is not constrained by the GERDA limit on either side. 
These results can easily be tested by GERDA in the near future. 




FIG. 11: The variation of Am§ 2 with respect 
to M ee for Case 2 of HSMR. 


FIG. 12: The variation of 023 with respect 
to 0i3 for Case 2 of HSMR. 


We show the variation of 023 with respect to 0 i 3 in Fig.[12j As can be seen, for the lowest 
possible value of a = 0.45, the allowed range is just a point which is located at 0 i 3 = 7.65° 
and 0 2 3 = 52.5°. As a = 1 . 1 , the range of 0i 3 is 7.65° — 8.4° and that of 0 2 3 is 48.5° — 52.5°. 
Finally for the highest value of a = 2. 5, 013 has almost the whole 3 a range 7.92° — 9.88° and 
the range of 0 2 3 is 36.8° — 48°. These results can be contrasted to case 1 where the minimum 
of the mixing angle 023 = 44.04° also happens for the upper value of a namely a = 1.28 and 
in both cases, the value of 023 can be in the first octant, contrary to the HSMU hypothesis. 
Also in both cases, the maximum of the mixing angle of 623 corresponds to the lower value 
of alpha. 

The next to be considered is the variation of 0 i 2 versus M ee as is shown in Fig. [T2J In 
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FIG. 13: The variation of $12 with respect 
to M ee for Case 2 of HSMR. 


FIG. 14: The variation of Em, with respect 
to M ee for Case 2 of HSMR. 


the right panel it can be observed that a = 0.45 corresponds to the minimum #12 = 30.8° 
and 0.384 eV < M ee < 0.42 eV. For a = 1.1 on the left panel, the whole 3 a range for 6*12 
is allowed for 0.39 eV < M ee < 0.4 eV and for M ee < 0.39 eV the allowed range of 2 
decreases. The upper value of a = 2.5, on the left panel, corresponds di 2 = 36.02° while 
0.344 eV < M ee < 0.366 eV. 

In Fig. [14l we show the behavior of the sum of the neutrino masses with respect to M ee . 
As can be seen in the right panel, for the lowest value of a = 0.45, Em,; lies in the range 
1.16 — 1.2 eV which corresponds to 0.383 eV < M ee < 0.4 eV. For M ee > 0.4, the range of 
Em, is 1.2 — 1.26 eV. For the upper value of a = 2.5, we have Em; = 1.05 — 1.12 eV while 
0.342 eV < M ee < 0.366 eV. On the left panel, for a = 1.1, the range of Em ; is 1.13 —1.2 eV 
which corresponds to 0.36 eV < M ee < 0.4 eV and the rest of the data point corresponds to 
M ee >0.4 eV. 


5. Case 3: 0\2 = 0\ 2 , $13 = 6 ^ 3 , $23 = « $23 

We now consider the final case where two of the leptonic mixing angles # 12 , $13 are 
identical to the quark mixing angle d q i2 , and the third leptonic mixing angle #23 is 
proportional to the quark mixing angle 0 23 . The correlation between A m 32 and M ee is 
shown in Fig. [13 The minimum allowed value of a, with all the mixing parameters within 
the global range, is 0.324. However in this case as can be seen from the right panel of 
Fig. E>1 we have 0.62 eV < M ee < 0.66 eV. This value of a corresponds to the entire 
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3a range of Am | 2 and violates the upper limit from GERDA. Therefore we also consider 
the lower value of a = 0.89, with all the mixing parameters within the global range and 
0.4 eV < M ee < 0.43 eV. This value corresponds to the entire 3a range of Am| 2 . The 
prediction of M ee at a = 1.1 is 0.372 eV < M ee < 0.41 eV corresponding to the entire 3a 
range of Am| 2 - The upper allowed value of a in this case is 1.52, (left panel of Fig. fT5l) 
with 0.3 eV < M ee < 0.34 eV and having the entire 3a range of Am§ 2 . Hence, the allowed 
range of a in this case covers the entire 3a range of A m\ 2 - The absolute neutrino mass scale 
increases for a < 1 and decreases for a > 1 similar to Case 2. The behavior of a in this 
case, Fig. [15] is different from Case 1, Fig. [TJ with the lower end of a constraining M ee . In 
this case it is possible to reach values of M ee as low as 0.3 eV compared to Cases 1, 2 and 
the HSMU hypothesis, and will only be ruled out if the limit from GERDA reaches M ee < 
0.3 eV. 
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FIG. 16: The variation of 0 2 3 with respect 
to 0i3 for Case 3 of HSMR. 

We next show the correlation of 0 2 3 with respect of $13 in Fig. [16] The 0 2 3 and $13 reach 
their 3a global limits at the lowest and upper most end of a. The value of 0 13 is 8.56° and 
that of 0 2 3 is 52.4° for lower allowed value of a = 0.89, corresponding to M ee = 0.4 eV. The 
allowed ranges of 0 i 3 and 0 2 3 for a = 1.1 in this case are much more constrained compared to 
Cases 1, 2 and the HSMU hypothesis. They are 7.62° — 8.05° and 50.1° — 52.1° respectively. 
The upper end of a = 1.52, results in a minimum value of 0 i 3 , whereas 0 2 3 is at maximum 
with 0 23 = 52.4°. The behavior of a here is different from Cases 1 and 2 with the lower end 
of a resulting in the upper end point of 0 i 3 and lower end point of 0 2 3 - 
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FIG. 15: The variation of Am| 2 with re¬ 


spect to M ee for Case 3 of HSMR. 
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FIG. 17: The variation of 9i 2 with respect 
to M ee for Case 3 of HSMR. 


FIG. 18: The variation of Em, with respect 
to M ee for Case 3 of HSMR. 


In Fig. [TT] we show the variation of 9 12 with M ee . We observe that the lower (30.78°) 
and upper (36.7°) 3cx global limits of Qyi , are reached at the upper most and the lowest ends 
of a respectively. In case of a = 0.89, the value of 0 12 is 30.62° which belongs to M ee = 0.4. 
The whole 3cr global range of 6* 12 is allowed for a = 1.1. 

Finally we show in Fig. [T8l the variation of the sum of the neutrino masses with respect 
to M ee . The region with M ee > 0.4 eV for a = 0.324, has Em, in the range of 1.84 — 2 eV 
and for a = 0.89 it is in the range 1.2 — 1.38 eV. In case of a = 1.1, with M ee < 0.4 eV, 
Em, is in the range 1.12 — 1.2 eV. The upper end of a = 1.52 has the sum in the range of 
0.93 — 1.01 eV. It is seen that Em; and M ee is much more relaxed compared to the HSMU 
and Cases 1, 2. 


6. Case 4 • #12 = a 0 q l2 , 9 13 = a 6^ 3 , 0 2 3 = 6t] 3 

We now consider the case where the leptonic mixing angles 6* 12 , 0 13 are proportional to 
the corresponding quark mixing angles 6 \ 2 , and the leptonic mixing angle 623 is identical 
to the quark mixing angle 0l 3 . The lowest allowed value of a for Case 4 is 0.92 which is 
derived using the 3cr global limits on mixing angles. The upper allowed value of a, respecting 
the GERDA limit, is 1.67. When we relax the GERDA limit then a turns out to be 1.77 
satisfying the 3cr global limits. We show, the correlation between A m \ 2 and M ee in Fig. 02B 
The lowest value of a = 0.92, covers the range 0.38 eV < M ee < 0.4 eV which corresponds 
to A m \ 2 = (2.30 — 2.50) x 10~ 3 eV 2 (cf. left panel of Fig. [19]). The prediction of M ee for 
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a = 1.1 is 0.384 eV < M ee < 0.41 eV corresponding to the whole 3a global range of Am 2 2 . 
On the right panel, the upper allowed end of a = 1.67 has 0.4 eV < M ee < 0.42 eV with 
Amjj 2 = (2.30 — 2.49) x 1CT 3 eV 2 . The upper most end a = 1.77 where the GERDA limit is 
not satisfied, has 0.46 eV < M ee < 0.48 eV corresponding to A mn\ 2 = (2.37—2.53) x 10 ~ 3 eV 2 , 
as shown in the right panel of Fig. [19j The behavior of a is similar to Case 1 with the upper 
values of a being constrained by the GERDA limit. The first distinction that this case 
offers, with the others considered before is that the whole 3cr range of Am 2 2 is not covered 
in Case 4 for all the allowed values of a. 



FIG. 19: The variation of Am| 2 with re¬ 
spect to M ee for Case 4 of HSMR. 


FIG. 20: The variation of 023 with respect 
to 0i3 for Case 4 of HSMR. 


Next, we show the correlation between 0 2 3 and $13 as illustrated in Fig. [201 The lower 
end of a = 0.92 reaches to the minimum of its 3a global limit for $13 and the maximum of 
the 3a limit for 0 2 3 - The situation for the upper most end of a = 1.77, is just opposite to 
the lower end, i.e. $13 is at the maximum of the 3a global limit whereas 0 2 3 is at its global 
minimum. This observation is just opposite to Case 3, where $13 ($ 23 ) reaches the global 
minimum (maximum), at the upper end of a. The allowed ranges of $13 and 0 23 , at a = 1.1 
for this case are 7.62° — 9.2° and 45.41° — 52.17°, respectively. The value of 0 13 is 9.59° and 
that of 0 23 is 37.71° - 37.76°, for a = 1.67. 

The variation of 0i 2 with respect to M ee is shown in Fig. [2U The lower and upper 3a 
global limits of 0i 2 are obtained at the lower and upper most end of a respectively. This 
observation is in contrast with Case 3 and Case 5 (to be discussed later). We get the full 
range of 0i 2 (cf. Fig. [21] for details) for a = 1.1, with some higher values of 0 12 being ruled 
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FIG. 21: The variation of #12 with respect FIG. 22: The variation of Em, with respect 
to M ee for Case 4 of HSMR. to M ee for Case 4 of HSMR. 

out for M ee < 0.394 eV. The value of 612 at a = 1.67 is 30.59° — 30.66°. 

Finally, we show the variation of the sum of neutrino masses with respect to M ee in 
Fig. [221 We find that for the lowest value of a = 0.92, the sum of neutrino mass ranges 
between 1.15 — 1.20 eV for 0.38 eV < M ee < 0.4 eV. In case of a = 1.1, Emj has a range 
of 1.157 — 1.23 eV corresponding to 0.38 eV < M ee < 0.41 eV as can be seen from the left 
panel of Fig. |22j A close look at the right panel of Fig. [22] reveals that for a = 1.67, Em, is 
in the range 1.22 — 1.27 eV corresponding to 0.4 eV < M ee < 0.42 eV. The range for Em* 
turns out to be 1.4 — 1.45 eV for a = 1.77 which corresponds to M ee > 0.4 eV. 


7. Case 5: 9 i2 = 0 q l2 , #13 = a 6 ^ 3 , 6*23 = a #| 3 

We now look at the case of the leptonic mixing angle 6* 12 being identical with its CKM 
counterpart and the other two leptonic mixing angles being proportional to the quark mixing 
angles. The correlation between Am 32 and M ee is shown in Fig. [221 as a deviates from unity. 
The minimum allowed value of a, with all the mixing parameters within the global range, 
is 0.06. However in this case as can be seen from the right panel of Fig. [23] we have 
2.24 eV < M ee < 2.28 eV, which violates the upper limit from GERDA. Therefore including 
the constraints of GERDA, the lowest possible value of a becomes 0.89. For a = 0.89, as can 
be seen from the right panel of Fig. [221 we obtain 0.40 eV < M ee < 0.42 eV corresponding 
to Am 32 = (2.20 — 2.48) x 10 ~ 3 eV 2 . The prediction of M ee for a=l.l from the left panel of 
Fig. E2]is 0.36 eV < M ee < 0.41 eV which belongs to the whole 3a range of Am| 2 - The upper 
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allowed value of a in this case is 3.18, (left panel of Fig. I23|) with 0.214 eV < M ee < 0.223 
eV corresponding to A m \ 2 — (2.28 — 2.46) x 10 ~ 3 eV 2 . The absolute neutrino mass scale 
increases for a < 1 and decreases for a > 1 similar to Cases 2 and 3. The behavior of a 
in this case is similar to Cases 2 and 3, with the lower end of a being constrained by the 
GERDA limit. It can also be seen from Fig. [23] that as we move towards the upper and 
lower ends of a, the whole 3a range of Am § 2 is not covered. 





FIG. 23: The variation of AmJ 2 with re- FIG. 24: The variation of 0 2 3 with respect 
spect to M ee for Case 5 of HSMR. to $13 for Case 5 of HSMR. 


We next show the correlation of 0 2 3 with respect of 0 13 in Fig. [24]. The 3a global end 
point limits of 623 and 0 13 , are reached at the lowest and upper ends of a. The values of 0 13 
and d 2 3 for lowest value of a = 0.89 are 8.40° and 52.23° respectively. These values belong 
to M ee = 0.4 eV. The allowed ranges of 0i 3 and 0 23 , at a = 1.1 for this case are 7.62° — 8.31° 
and 49.0° — 52.3° respectively. The upper end of a = 3.18, results in a global minimum 
value of 613 and a global maximal value of 0 23 , similar to Case 3. The lower end of a results 
in a global maximum value of 0 13 and a global minimum value of # 23 . 

We next show the variation of 0 12 with M ee in Fig. [25J The lower (30.62°) and the upper 
(36.81°) 3a global limits of # 12 , correspond to the upper and the lowest ends of a. However, 
for case of a = 1.1, the whole 3a range of 0\2 (30.62° — 36.81°) is covered. 

Finally we show in Fig. [26] the variation of the sum of the neutrino masses with respect 
to M ee . The region with M ee > 0.4 eV, for a = 0.06, has Em* in the range of 6.73 — 6.85 eV 
and for a = 0.89 it is in the range 1.20 —1.28 eV. For a = 1.1, with M ee < 0.4 eV, Em* is in 
the range 1.07 — 1.22 eV. The upper end of a = 3.18 has the sum in the range of 0.65 — 0.67 
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FIG. 25: The variation of #12 with respect 
to M ee for Case 5 of HSMR. 


FIG. 26: The variation of Smj with respect 
to M ee for Case 5 of HSMR. 


eV which is remarkably lower than the previous cases. We note that the Em, is below the 


cosmological upper bound 


52], The further discussion on the cosmological constraints on 


our work will be provided in the last section of this paper. It is observed that Case 5 behaves 
almost similarly to Case 3, and is the most relaxed one in terms of M ee . We can go to values 
of M ee as low as 0.21 eV, consistent with the upper end of a. Hence, this case is partially 
beyond the reach of GERDA sensitivity which is maximum 0.3 eV. However, this is well 
within the reach of KATERIN experiment [dS ], 


8 . Case 6: #12 = a 0\ 2 , $13 = # 13 , #23 = a #23 

We next consider the case where the leptonic mixing angle $13 is identical with its CKM 
counterpart and the other two leptonic mixing angles are proportional to the quark mixing 
angles. The correlation between A m \ 2 and M ee is shown in Fig. [23 The minimum allowed 
value of cc, with all the mixing parameters within the global range, is 0.86. It can be seen 
from the left panel of Fig. [23 for this value of a, M ee has a range 0.397 eV < M ee < 0.42 eV 
corresponding to A m \ 2 = (2.20 — 2.48) x 10~ 3 eV 2 . In case of a = 1.1, we have 0.36 eV < 
M ee < 0.424 eV which corresponds to the whole 3cr range of A m\ 2 . At the upper allowed 
value of a = 2.11 as seen from the right panel of Fig. [23 we have 0.41 eV < M ee < 0.45 eV 
with the whole 3 a range of A m \ 2 covered. We have 0.64 eV < M ee < 0.67 eV corresponding 
to A 77132 = (2-22 — 2.54) x 10” 3 eV 2 for the uppermost value of a = 2.19. This end is 
already rejected by the GERDA limit. In this case, it is worth mentioning that the absolute 
















29 


neutrino mass scale increases for both a < 1 and a > 1. 




FIG. 27: The variation of Amj 2 with re- FIG. 28: The variation of #23 with respect 
spect to M ee for Case 6 of HSMR. to #13 for Case 6 of HSMR. 


The behavior of a in this case Fig. [23 is different from Case 5, Fig. [221 Unlike Case 5 
considered before, in this case for both lower and upper end values of a we get M ee close 
to its upper limit. This is because unlike the previous cases, in this case the limit on lower 
value of a comes not from M ee but from the neutrino oscillation parameters. Also, although 
the upper limit of a = 2.11 is constrained by M ee but this value is quite close to the upper 
limit of 2.19 obtained without the M ee constraint. 

In Fig. [221 we show the correlation of 6*23 with respect of # 13 . The 3a global end point 
limits of $23 and # 13 , are reached at the lowest and uppermost end of a. The allowed ranges 
of #13 and # 23 , at a = 1.1 for this case are 7.62° — 8.90° and 47.0° — 52.3° respectively. At 
a = 2.11, the value of #13 is 9.4° — 9.76° and that of #23 is 38.7° — 38.8°. 

We next show the variation of #12 with M ee in Fig. [2HJ The lower (30.60°) and upper 
(36.81°) 3a global limits of # 12 , is reached at the lowest and the uppermost end of a. This 
behavior is quite the opposite of the behavior shown in Fig. [251 for Case 5. In case of a = 
1.1, the whole 3a range of #i 2 (30.60° — 36.81°) is covered. The value of #i 2 at a = 2.11 is 
34.57° - 35.02°. 

Finally we show in Fig. [301 the variation of the sum of the neutrino masses with respect 
to M ee . I 11 case of a = 0.86 the sum of neutrino masses Em; is in the range of 1.19 — 1.27 
eV. Next for a = 1.1, Ewj has a range of 1.1 — 1.27 eV for M ee < 0.4 eV and when a = 2.11 
it is in the range 1.26 — 1.36 eV for 0.4 eV < M ee < 0.48 eV. 
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M ee < eV > 





FIG. 29: The variation of #12 with respect FIG. 30: The variation of Em, with respect 
to M ee for Case 6 of HSMR. to M ee for Case 6 of HSMR. 

9. Case 1: 0\2 = a 0\ 2 , # 1,3 = a O 23 = a #23 

We finally consider the case where all the leptonic mixing angles are proportional to the 
quark mixing angle by the same proportionality constant ( a ). We find that the upper bound 
on a, is constrained by the mass limit (M ee ) from GERDA, whereas the lower limit on a is 
constrained by the 3a global limit of the leptonic mixing angles. The lowest value of a is 
0.89 and the highest value of a relaxing the GERDA limit is 2.09, whereas by taking into 
account the M ee limit, the highest value is 2. 

We next discuss the behavior of the neutrino mass and mixing parameters in Case 7, with 
the variation of a in the allowed range. Firstly like all the previous cases, the variation of 
Ar?7| 2 with M ee is shown in Fig. EU As seen from the left panel of Fig. [3TJ for the lowest 
value of a = 0.895, we have 0.391 eV < M ee < 0.425 eV which corresponds to the whole 
3cr global range of Am§ 2 . At a = 1.1, the range of M ee is 0.362 eV < M ee < 0.405 eV 
which again corresponds to the whole 3a global range of A m\ 2 . The range of M ee at a = 2 
is 0.4 eV < M ee < 0.42 eV which corresponds to the whole 3cr global range of Am| 2 . The 
uppermost end of a = 2.09 has 0.42 eV < M ee < 0.452 eV which corresponds to the whole 
3cr global range of A m \ 2 and is rejected by the GERDA limit. Hence, the entire allowed 
range of a covers the whole 3cr range of m\ 2 - The absolute neutrino mass scale increases for 
both a < 1 and a > 1 similar to Case 6. The behavior of a resembles to Case 6 with the 
upper and lower ends of a having values close to M ee . 

We show the range of 623 and Q \ 3 covered by the different allowed values of a in Fig. [221 
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FIG. 31: The variation of Am| 2 with re¬ 
spect to M ee for Case 7 of HSMR. 


FIG. 32: The variation of 023 with respect 
to 0i3 for Case 7 of HSMR. 


The 3a global limits on the mixing angles are reached at the lower and uppermost ends 
of a. The allowed ranges of $13 and 023 for a = 1.1 are 7.76° — 9.02° and 46.3° — 52.26° 
respectively. For a = 2.0, the value of 0i 3 is 9.44° and that of 0 2 3 is 37.8° which belongs to 
M ee = 0.4 eV. The uppermost end of a = 2.09, gives the value of 0 13 at its global upper 
limit, whereas 023 is kept at its global lower limit. The converse is true for the lower end of 
a with 0 i 3 , 0 i 2 at its lower value and 023 at its maximum. 

The variation of the third mixing angle 0 i 2 with respect to M ee is next plotted in Fig. [33l 
The pattern obtained is similar to Case 6 , with the lower and upper end of a giving the 3a 
global end points of 0i 2 respectively. The whole 3a global range of 0 12 is allowed, for a = 
1.1. The value of 0 i 2 at a = 2 is 36.6°. 

Finally we plot the sum of the neutrino masses as a function of M ee in Fig. [34j For a = 
0.895, range of Sm,; is 1.18 — 1.28 eV corresponding to 0.391 < M ee < 0.425 eV. The range 
of Em* at a = 1.1 is 1.09 — 1.18 eV for M ee < 0.4 eV. At the upper allowed value of a = 2.0, 
it is 1.2 — 1.28 eV for 0.4 < M ee < 0.425 eV. The sum of the neutrino masses is 1.27 — 1.38 
eV for 0.42 < M ee < 0.452 eV, in case of the uppermost value of a = 2.09. This range is 
not allowed by the GERDA limit. 

Lastly as a completion, in order to give a clear picture of all the cases discussed here 
along with their phenomenological consequences, we summarize our results in Table 1III1 
The upper and lower ends of a allowed by the experiments for all the cases are presented 
along with the corresponding values of masses and mixing angles of the neutrino sector. 
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FIG. 33: The variation of #12 with respect FIG. 34: The variation of Smj with respect 
to M ee for Case 7 of HSMR. to M ee for Case 7 of HSMR. 

10. The effects of new physics within type-1 seesaw framework 

In this sub-section, we discuss the possible effects of the new physics which could generate 
dimensional-5 operator. For sake of illustration we take type-I seesaw as the mechanism 
responsible for generating the effective dimensional-5 operator. The RG equations for type-I 


seesaw can be found in the Ref. 53]. In table IIV1 we show results for different cases of 
HSMR with a seesaw scale equals 4 x 10 13 GeV which is slightly lower than the scale of 
the dimensional-5 operator. We have chosen this scale to demonstrate and differentiate the 
effects of type-I mechanism from the results obtained using dimesional-5 operator. Since 
the major part of RG magnification happens at scales much lower than the typical seesaw 
scales, the results obtained from dimension-5 operator and those obtained form type-I seesaw 
mechanism are not very different over a large range of parameters 


28 


32 


35|. Our results 


for type-I seesaw are as shown in Table [IVl To take the effect of type-I seesaw thoroughly, 
we have done the RG running from the GUT scale (2 x 10 16 ) GeV to the seesaw scale using 
the full RG equation for type-I seesaw mechanism. Below the seesaw scale the right handed 
neutrinos are integrated out and as before, the subsequent RG running is done with effect 
dimension-5 operator. Here we will like to remark that since in this case the RG running 
is done from a higher scale i.e. GUT scale so we expect small deviations from the previous 
results primarily due to the larger range of RG running. The dependence of RG evolution 


on the chosen high scale is studied in [32|, |35]. 


We observe from comparing Tables IIVI and IIIII that for case 1 of the HSMR, the lower 
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a 

Masses at unification scale (eV) 

S rm (eV) 

^12 

^13 

^23 

Am| a (10 -3 eV a ) 

Ami, (lCT 5 eV 2 ) 

M ee (eV) 

Lightest neutrino mass: mi (eV) 

R 

mi 

m 2 

m 3 

HSMU hypothesis 

1 

0.45700-0.47686 

0.46-0.48 

0.51573-0.53817 

1.16-1.2 

30.59-36.81 

7.63-8.34 

49-52.3 

2.21-2.45 

6.99-8.18 

0.384-0.4 

0.38397 

1.5-1.8 

Case 1 

0.902 

0.44214-0.47791 

0.445-0.481 

0.49891- 0.53927 

1.12-1.2 

30.6 

7.62 

52.29 

2.20-2.55 

6.99-8.18 

0.365-0.4 

0.40155 

1.39-1.42 

1.28 

0.47573 

0.47849 

0.53772 

1.2 

32.82 

8.16 

44.05 

2.20 

7.196 

0.4 

0.39968 

2.1 

Case 2 

0.45 

0.41 - 0.445 

0.411 - 0.448 

0.464 - 0.54 

1.16-1.2 

30.8 

7.65 

52.3 

2.20-2.57 

7.00-8.18 

0.38-0.42 

0.345-0.375 

1.70-1.95 

2.5 

0.415 - 0.444 

0.416 - 0.445 

0.468 - 0.5 

1.05-1.12 

31.89-36.6 

7.92-9.88 

37.7-48 

2.20-2.53 

7.00-8.18 

0.342-0.378 

0.348-0.373 

1 

Case 3 

0.89 

0.4771 

0.4800 

0.5384 

1.2 

30.6 

8.56 

52.3 

2.20 

6.99 

0.4 

0.4009 

1.44 

1.52 

0.3693-0.3990 

0.3730-0.4030 

0.4179-0.4515 

0.93-1.01 

30.78 

7.62 

52.3 

2.20-2.56 

6.99-8.18 

0.31-0.336 

0.3102-0.3352 

2.16-2.21 

Case 4 

0.92 

0.4550-0.4749 

0.458-0.478 

0.5135-0.5359 

1.15-1.20 

30.6 

7.65 

52.28 

2.30-2.50 

6.99-8.01 

0.382-0.40 

0.3823-0.3990 

1.463-1.479 

1.67 

0.4819-0.5018 

0.485-0.505 

0.5460-0.5685 

1.22 -1.27 

30.59-30.66 

9.59 

37.71-37.76 

2.29-2.49 

7.08-8.18 

0.405-0.422 

0.4049-0.4216 

2.82-2.86 

Case 5 

0.89 

0.475-0.505 

0.478-0.508 

0.536-0.570 

1.20-1.28 

30.72 

8.40 

52.23 

2.20-2.48 

7.00-7.97 

0.4-0.42 

0.399-0.424 

1.465-1.485 

3.18 

0.255-0.265 

0.26-0.27 

0.29-0.30 

0.65-0.67 

30.62 

7.62 

52.3 

2.27-2.45 

6.99-8.12 

0.214-0.223 

0.214-0.222 

3.542-3.672 

Case 6 

0.86 

0.47-0.50 

0.474-0.504 

0.531-0.565 

1.19-1.27 

30.61 

7.62 

52.3 

2.20-2.48 

6.99-8.17 

0.396-0.421 

0.3961-0.4211 

1.241-1.268 

2.11 

0.494-0.533 

0.502-0.542 

0.567-0.612 

1.26-1.36 

34.57-35.02 

9.64-9.76 

37.70-37.8 

2.20-2.57 

6.99-8.18 

0.416-0.450 

0.414-0.448 

246.942-251.842 

Case 7 

0.895 

0.4683-0.5051 

0.4710-0.5080 

0.5282-0.5696 

1.185-1.22 

30.68 

7.63 

52.28 

2.20-2.31 

6.99-8.18 

0.393-0.4 

0.3935-0.4044 

1.34-1.37 

2 

0.4738 

0.4800 

0.5413 

1.18 

36.6 

9.44 

37.8 

2.20 

7.248 

0.39 

0.3977 

81.7-82.11 


TABLE III: The allowed predictions for HSMU and the different cases of the HSMR for lower and 
upper allowed values of a, Eqs. (I14II20I) . 

allowed end of a effectively does not change. As expected, there are slight changes in the 
value of the observables. For example, the M ee decreases and reaches to the value 0.349 
eV compare to the prediction given in table IIIII Similar observation for the mass of the 
lightest neutrino. The upper end of a changes after introducing seesaw scale, primarily due 
to increased RG running range. In Table [Till the upper allowed end for case 1 is 1.28. As 
we observe in Table El it is now 1.71 and parameter space is bit expanded. However, there 
is no significant qualitative change in our results which are same as before. 

Similarly, for case 2, one can observe from the Tables ITVl and m that the lower end of a 
does not change much. The results are stable and similar to Table IIIII The upper end of a 
changes slightly and is 2.59 now. Again as before, there is no significant qualitative change 
in our results which are same as before. 

In case 3 of Table IIV1 the lower value of a has shifted a bit from that of Tabl dUIl that 
is from 0.89 to 0.75, but the higher value remains intact being 1.52. The parameters also 

































34 


cover more or less the same span as before and show a stable situation . 

The observed pattern in case 4 is same as the case 1, as we see that the lower end does not 
change while the upper end changes from 1.67 to 1.76 after the inclusion of type-I seesaw 
(see Table hyi and Table nn for comparison). The value of M ee and the mass of lightest 
neutrino decreases compare to the values given in Table Hill and attain the values 0.3508 eV 
and 0.3506 eV, respectively. 

The observed pattern in case 5 is also same as that obtained in Table IIIII Comparing the 
results with Table m we find that the value of a at the lower end changes slightly. The 
lower end saturates the bound for M ee , whereas at the upper end the value of M ee turns out 
to be 0.205 eV, which is slightly smaller than the value quoted in Table HTTl 

In case 6, the upper value of a changes very slightly compared to Table IIIII whereas 
the lower value remains the same. In this case both the lower and upper end saturates the 
bound for M ee . 

As expected, the results for case 7 are stable as can be observed from Tables Hill and HVl 
where the lower end of a is now 0.896 instead of 0.895 and the higher end remains the same, 
namely a = 2. 

Thus, as expected the results obtained with the framework of type-I seesaw mechanism 
are qualitatively same as those obtained using only dimension-5 effective operator. The 
general observation is that the absolute mass scale is decreasing due to the RG evolution 
starting from GUT scale (2 x 10 16 GeV) which is higher than the scale of dimensional-5 
operator. This leads to a slight change in the allowed end of a that is constrained by the 
observable M ee . That is why, we observe a slight change in the values of mixing angles. We 
remark that if the high scale from where RG evolution begins, is chosen to be 10 14 GeV 
with a seesaw scale equals 4 x 10 13 GeV, we recover the results obtained with dimensional-5 
operator which is naturally expected. 


V. THEORETICAL MODELS FOR HIGH SCALE MIXING RELATIONS 


In this section, we address the theoretical implementation of HSMR hypothesis from the 
model building point of view. The only aim of this section is to illustrate that the HMSR 
hypothesis can be simply realized in models based on flavor symmetries. We follow the same 
line of argument as presented in Ref. 
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a 

Masses at unification scale (eV) 

Emi (eV) 

e \ 2 

e\ 3 

023 

Amj 2 (10 3 eV 2 ) 

A(10 _ 5 eV 2 ) 

M ee (eV) 

Lightest neutrino mass: mi (eV) 

R 

mi 

m 2 

m 3 

°“1 

0.903 

0.46547 

0.469 

0.53833 

1.04 

31.06 

7.63 

52.30 

2.209 

8.00 

0.349 

0.34419 

1.48 

1.71 

0.54038-0.54199 

0.545-0.546 

0.62793-0.62968 

1.203- 1.205 

30.82 - 33.38 

8.68 - 9.87 

37.75 - 40.25 

2.201-2.26 

7.05-7.96 

0.399-0.4 

0.39815 

4.89 

Case 2 

0.47 

0.486067-0.486077 

0.49 

0.56241 - 0.562434 

1.08-1.0821 

30.66-30.92 

7.62 

52.15-52.19 

2.2-2.23 

7.65-8 

0.359 

0.359022-0.359308 

1.96-1.99 

2.59 

0.478-0.480 

0.481- 0.483 

0.55 - 0.554 

1.06 

30.68-31.2 

9.1-9.16 

42.96-43.51 

2.5-2.56 

7-7.24 

0.354 

0.353-0.354 

1 

Case 3 

0.75 

0.541797 - 0.542934 

0.545 - 0.546 

0.625782 - 0.627092 

1.20171 - 1.20407 

30.6074 - 32.7328 

8.74619 - 8.95191 

51.3341 - 52.2853 

2.2-2.24 

6.99-8.16 

0.399 - 0.4 

0.399 - 0.4 

1.351 - 1.374 

1.52 

0.392477 

0.3971 

0.455448 

0.879869 

30.7133 

7.63895 

52.2976 

2.28 

7.05 

0.291584 

0.291263 

2.52 

°” 4 

0.92 

0.4744 

0.478 

0.5487 

1.06 

30.98 

7.64 

52.27 

2.266 

8.081 

0.3508 

0.3506 

1.55 

1.76 

0.5484-0.5544 

0.552-0.558 

0.6366-0.6435 

1.22-1.23 

36.09-36.55 

9.80-9.87 

37.71-37.83 

2.204-2.207 

7.101-7.837 

0.4045-0.4088 

0.4039-0.4083 

4.7-4.85 

C"”5 

0.80 

0.537 

0.540 

0.621 

1.191 

32.68 

7.79 

49.17 

2.20 

8.09 

0.396 

0.396 

1.532 

3.18 

0.274 

0.280 

0.321 

0.623 

30.65 

7.63 

52.3 

2.45 

7.22 

0.205 

0.205 

4.582 

Case 6 

0.86 

0.497 

0.500 

0.574 

1.104 

31.09 

7.62 

52.3 

2.20 

7.86 

0.367 

0.367 

1.301 

2.14 

0.544 

0.555 

0.642 

1.22 

35.06 

9.88 

37.77 

2.21 

8.18 

0.402 

0.400 

957 

Case 7 

0.896 

0.49 

0.50 

0.57 

1.104 

30.65 - 31.41 

7.62 

52.26 

2.24 

7.1-8.1 

0.36 

0.36 

1.42 

2 

0.5 

0.49-0.51 

0.58 

1.12 

34.30 - 36.69 

9.52 - 9.84 

37.80 - 38.07 

2.28 - 2.47 

7.08 - 7.21 

0.37 

0.37 

155.5 


TABLE IV: The allowed predictions for the different cases of the HSMR for lower and upper 
allowed values of a, Eqs. (|14H2(JI) within the framework of type-I seesaw for sea-saw scale 4 x 10 13 
GeV. It should be noted that RG evolution begins from GUT scale which is 2 x 10 16 GeV. 

Now we discuss a simple realization of HSMR relations using abelian Z 7 flavor symmetry. 
To realize the HSMR relations we add three SU (2) triplet scalars i = 1, 2, 3 to the particle 
content of MSSM. The smallness of neutrino masses can then be explained by the type-II 
seesaw mechanism. Let the quarks and leptons and scalars transform under Z 7 as follows 

Q\ ~ 1) Ql ~ u i Ql ~ a;3 ) u R ,d R ~ 1, c R ,s R ~ uj, t R ,b R ~ u 3 
L\ ~ 1, L\ ~ u, L{ ~ cu 3 , e R ~ 1, n R ~ u, r R ~ a; 3 

H u ,H d ~ 1 , ~ 1 , £ 2 ~ a ; 2 , £3 ~ uj 6 ( 46 ) 

where uj = e^ 1 is the seventh root of unity. In the above equation Q l L , U L ; i = 1, 2, 3 are 
the quark and the lepton doublets respectively whereas u R , d R , c R , s R , t R , b R , e R , /i R , t r are 
the quark and the charged lepton singlets. Moreover, H Ul H d are the two scalar doublets 
required to give mass to the up and down type quarks respectively. 

It is easy to see from Eq. (HR]) that the Z 7 symmetry leads to diagonal mass matrices for 
both the quarks and the leptons leading to Uckm = Upmns = I■ To obtain the realistic 
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CKM and PMNS matrices as well as the HSMR relations, we allow for small Z 7 symmetry 


breaking terms as done in Ref. 


541 ] albeit for A 4 symmetry. Such corrections can arise 


from soft supersymmetry breaking sector as shown in Ref. 


55m71]. Allowing for symmetry 


breaking terms of the form | h'"\ « \h"\ « \h[\ « \hi\ where hi are the terms invariant 
under Z 7 symmetry and h![ and h"' are the symmetry breaking terms transforming as 
oj, u} 2 and a; 3 respectively under Z 7 symmetry. Following the approach of Ref. 54] one can 
then easily realize the HSMR relations. Here we want to emphasize that owing to quite 
different masses of quarks and charged leptons, this analysis will in general lead to HSMR 
relations and not to HSMU relations. To obtain HSMU relations from such an approach 
one has to invoke a symmetry or mechanism to ensure that the symmetry breaking terms 
are exactly same in both quark and lepton sectors. 

Before, ending this section we would like to further remark that although we have only 
discussed realization of HSMR relations through the Z 7 symmetry, they can also be quite 
easily and naturally realized using other flavor symmetries and also using other type of 
seesaw mechanisms. For example, one can also realize HSMR relations within the framework 
of type-I seesaw mechanism using Z 7 symmetry. For this, instead of adding triplet scalars 
we add three right handed neutrinos which transform as N 4 r ~ 1, N 2 r ~ oj, N 3 r ~ ca 3 under 
the Z 7 symmetry. We also add three heavy singlet scalars 0*; i = 1,2,3 transforming as 
4>i ~ 1, 02 ~ co 5 , 03 ~ cu under Z 7 symmetry. Following computations analogous to those 
done above, one can again easily obtain the HSMR relations. Thus, it is clear that HSMR 
relations are very natural and can be easily realized using discreet flavor symmetries. In 
this work we do a model independent analysis of the consequences of the HSMR relations 
assuming they are realized at the high scale by appropriate flavor symmetries. 


VI. SUMMARY 

The very small mass of the neutrinos along with a large mixing among them is arguably 
a remarkable observation. This phenomenon is starkly different from the mixing in the 
quark sector which is small in the SM. The quest to understand the origin of a large mixing 
among the neutrinos and a small mixing among the quarks has led to many interesting 
theoretical ideas. Many beyond the standard model scenarios have been constructed, trying 
to understand the major theoretical challenge posed by the neutrino mixing. GUT theories 
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with the quark-lepton unification have been extensively used in the literature to understand 
the neutrino sector at low energies. The postulated HSMR in this work is another effort 
to understand this extraordinary observation of neutrino mixing. We have shown from a 
model building point of view, how the HSMR can be naturally realized using different flavor 
symmetries and seesaw mechanisms. We have first considered the most general relation 
among the leptonic and the quark mixing angles, with different proportionality constants 
(cq). We then list the different possible cases which arise, for the maximum and minimum 
allowed values of a*. It is found that for the allowed range of ctj, M ee is between 0.35 eV 
- 0.4 eV. The future experiments from GERDA will severely constrain these scenarios. We 
then look into more simplified cases to have a clear physical picture and therefore consider 
the cp to be equal for the three generations and vary A;* to 0 or 1. We then list the seven 
possible ways the quark and the leptonic mixing angles can be proportional to each other (cf. 
Eqs. ([T41 - I2U1A It is remarkable that these relations naturally explain the difference between 
Vckm and Upmns at the low scale. Furthermore, the QLC relation and the observation in 
Eq. ([2]) can be easily recovered by these relations. 

We have thoroughly investigated the implications and the phenomenological consequences 
of all the possible cases, taking into account the latest experimental constraints. The whole 
analysis has been done with the assumption of normal hierarchy and QD mass pattern. 
In general, we have discovered three new correlations among A m\ 2) M ee , 0 12 and sum of 
neutrino masses. These correlations are not investigated in previous studies. 

We first discuss about the HSMU scenario, which is a special case of all the HSMR 
scenarios in the a = 1 limit. The behavior of the neutrino masses and the mixing parameters 
at the low energy scale is discussed in detail for all the cases in HSMR with the value of a 
deviating from unity in the allowed range. The allowed range of a is bounded by the recent 
experimental results listed in Table Q] and the upper limit on M ee provided by GERDA 511. 
It is seen that for all the cases except Case 2, the M ee constraint from the GERDA results in 
either upper (Cases 1, 4, 6 and 7) or lower (Cases 3 and 5) limit of a. Otherwise the allowed 
value of a is mostly constrained by the 3<r global limits on neutrino mixing parameters. 

An interesting feature is observed in Case 2, where the lower end is constrained by the 3 a 
global limits of neutrino mixing parameters but the upper end is constrained by the value of 
the ratio R, which contributes through threshold corrections. We have worked here in the 
inverted hierarchy scheme in the charged slepton sector, forcing the ratio to be either greater 
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than or equal to one. A common behavior has been observed for all the cases, where we 
always find a strong correlation between 6*23 and # 13 , for all the allowed values of a except at 
the end points which corresponds to a point in the 623 — 613 plane. It is also seen that among 
all the experimental constraints M ee is the most interesting one as it mostly constrains the 
different cases as well as differentiates among them. If in the future the upper limit from 
GERDA goes down to 0.35 eV, then HSMU, Case 1, Case 4, Case 6 and Case 7 will be 
ruled out. The ones who will survive will be Cases 2, 3 and 5 which allows M ee as low as 
0.2 eV but with the value of a > 1. The constraint on M ee can automatically be reverted 
to the sum of the neutrino masses. It will show a similar behavior while discriminating the 
various cases. We also notice that if we take into account the GERDA limit of 0.4 eV, then 
the allowed range of a in Cases 1, 3 and 5 is limited to a small region in the 623 — 613 plane 
(Figs. E DU [231)• Therefore these cases along with HSMU will be ruled out, if the best fit 
value of 623 becomes less than 44° or that of #13 becomes greater than 8.55° in the future. 
We further see that Cases 2 and 5 can survive longer, and the region of M ee = 0.2 eV is 
beyond the sensitivity of GERDA which is maximun 0.3 eV. The region of M ee = 0.2 eV will 


easily be probed by KATRIN 


48] since mp is approximately identical to M ee in this work. 


Here, we pause to comment on the cosmological limit on the sum of the neutrino masses 


52]. 


Our predictions in all cases except Case 5 are slightly above the upper cosmological bound 
of 0.72 eV. As commented earlier, this bound is model dependent. Hence, it is preferred to 
test predictions of this work in a laboratory based experiment, like GERDA 511] . 

We also observe that Cases 3 and 5 show similar behavior, this is mainly because both 
consider the framework, where the neutrino mixing angle 623 is equal to the quark mix¬ 
ing angle by a proportionality constant cx.6\ 3 . Although Case 5 also has the condition of 
6 13 = a d\ 3 , but at the GUT scale df 3 <C 6r, 3 , therefore the effect of 6\ 3 dominates. The 
same pattern can be observed for Case 1 and Case 4, explained through the same argument, 
6 \ 3 < &l 2 at the GUT scale. Continuing the same argument as expected we find that Case 
7 displays similar behavior as Case 6 . The effect of the neutrino mixing angle $13 being pro¬ 
portional to the quark mixing has many interesting results, as it leads to the most optimistic 
case. However once the other angles become proportional, this effect is subdued. Finally we 
note that all these interpretations have been done with the assumption that the Dirac and 
the Majorana phases of the PMNS matrix are zero and phenomenological consequences can 
change with nonzero phases. The overall scenario depicting a quark-lepton symmetry at a 
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high scale through HSMR can be narrowed down to a particular case or completely ruled 
out, only from the future improved experimental constraints. These constraints can be from 


the neutrinoless double beta decay 


511 ]. or the LHC constraints on the SUSY spectrum. 


The different scenarios of the HSMR can be discriminated through measurement of var¬ 
ious observables like M ee and by precise determination of the values of the mixing angles, 
particularly # 13 and 623 mixing angles. As we have shown in the figures for various cases 
as well in the tables, the allowed ranges for M ee and the angles are different for different 
cases and a precise determination of these observables can be used as a way to distinguish 
various cases of HSMR. In addition to neutrino observables one can also use other process 
like lepton-flavor violation to distinguish the different allowed cases. The mass-splitting in 
the charged-slepton sector is given by the ratio R = -rrfz. We observe from tables UT1 and 
mi that the ratio R almost discriminate every scenario. Hence the processes like p —>■ ey, 
/i —> eee and anomalous magnetic moment of the electron. For example, the SUSY contri¬ 
bution to the anomalous magnetic moment of the electron directly depends on the ratio R 


58l] . The detail study of this aspect of the work is not possible in this paper. 


Furthermore, for sake of completion, we also present our results in the framework of the 
type-I seesaw. The aim is to show how the predictions do not change in any significant way 
and that the analysis done with effective dimension-5 operator is quite robust. As argued 
before, this is not surprising as the major part of RG magnification happens only at much 
lower scales closer to SUSY breaking scale. At such low scales, the effective dimension-5 
operator provides a very good approximation to the high scale seesaw mechanisms. The 
mass scale of the right-handed neutrinos is chosen 4 x 10 13 GeV which is close to the scale 
of the dimensional-5 operator. We notice that parameter spcae increases very slowly as we 
decrease the scale of new physics primarily due to increased span of RG running. However 
predictions do not change in any significant manner and are quite robust. 

We also comment on a general theoretical view which is more general than the HSMU 
hypothesis and the HSMR. Assuming that at some high scale, both the mixing matrices 
(CKM and PMNS) are approximately unit matrices, but some perturbation can mix the 
generations leading to the Wolfenstein form of the mixing in both the quark and lepton 
sectors. This results in the mixing between the first and the second generations to be A (a 
small number of order 0.2), the second and the third generations mixing to be second order 
in A i.e. sin 623 ~ A 2 while the first and the third generations mixing to be third order order 
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in A i.e. sin 0 13 ~ A 3 . Now after RG evolution the CKM mixing angles do not change much 
but the PMNS mixing angles are dramatically magnified for the reasons already mentioned 
in this as well as our earlier papers 
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Finally, in short, crux of our paper is following. 

• We have proposed and studied the ffSMR hypothesis which is a more general frame¬ 
work than the ffSMU hypothesis. 


• The HSMR hypothesis provides a very simple explanation of the observed large neu¬ 
trino mixing. The present and future neutrino experiments can easily test predictions 
of our work. If our predictions are confirmed by experiments, like GERDA, it would 
be a good hint of quark-lepton unification at high scale. 


• We observe that the ffSMU hypothesis represents q: = 1 limit of the ffSMR hypothesis 
and is constrained by the lowest allowed value of M ee which is 0.384 MeV. Therefore, 
if the ffSMU hypothesis is ruled out by experiments, like GERDA, the other ffSMR 
cases with a ^ 1 may survive and their confirmation would be itself a strong hint of 
the proportinality between quark and leptonic mixing angles which is the basis of the 
ffSMR hypothesis. 


• We have done a rigorous, thorough and comprehensive study with the ffSMR hy¬ 
pothesis which does not exist in the literature. All results reported in the literature 
using the ffSMU hypothesis, are very small subset of our results with the ffSMR hy¬ 
pothesis presented in our paper. Moreover, we have also thoroughly compared ffSMR 
hypothesis with respect to the ffSMU. 


• In our work, we have discovered new strong correlations among different experimental 
observables for every limit of the HSMR hypthesis. These correlations do not exist 
in the literature and are easily testable in present ongoing experiments. For example, 
there is a strong correlation between A m\ 2 and M ee . This correlation can be easily 
tested by GERDA experiment. There are two more such correlations namely among 
9 12 , M ee , Y,rrii and M ee discussed in our work which are completely new and unexplored 
in the literature. 


• Furthermore, we have comprehensively studied a strong correlation between 0 2 3 and #13 
and predictions can be easily tested in present ongoing experiments. This correlation 
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was studied in a previous study in a specific limit. Since we have done a comprehensive 
full parameter scan, this correlation has become a robust band now. 
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